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Abstract 

Thermal and electric conductivities of relativistic degenerate electrons are calcu- 
lated for the case when electrons scatter by phonons in Coulomb crystals made of 
spherical finite-size nuclei at densities 10^^ g/crc? < p < 10^^ g/cm"^, corresponding 
to the inner crust of a neutron star. In combination with the results of the previous 
article (for lower p), simple unified fits are obtained which describe the kinetic coef- 
ficients in the range 10'^ g/cm'^ ^ P ^ 10^^ g/cm^, for matter with arbitrary nuclear 
composition. The results are valid for studying thermal evolution of neutron stars 
and evolution of their magnetic fields. The difference between the kinetic coefficients 
in the neutron star crust composed of ground state and accreted matters is analyzed. 
Thermal drift of the magnetic field in the neutron star crust is discussed. 
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1 INTRODUCTION 

The thermal end electrical conductivities of degenerate electrons which suffer electron- 
phonon scattering in Coulomb crystals of atomic nuclei in dense matter of white dwarfs 
and neutron stars has been studied by a number of authors. In particular, Yakovlev and 
Urpin (1980) developed an approximate method of evaluation of the kinetic coefficients 
which takes into account the main features of phonon spectrum in the body-centered-cubic 
Coulomb crystals and the leading role of the Umklapp processes in electron-phonon scat- 
tering under astrophysical conditions. The same authors presented critical analysis of the 
preceding papers. Raikh and Yakovlev (1982) confirmed the validity of the approximate 
method by direct Monte Carlo calculations of the kinetic coefficients with allowance for 
the exact phonon spectrum and contributions from the normal and Umklapp scattering 
processes. Itoh et al. (1984, 1993) used the same approximate method but included the 
Debye - Waller factor. These authors fitted their results by very complicated expressions. 
Finally, Baiko and Yakovlev (1995) calculated the kinetic coefficients directly by the Monte 
Carlo method and by the approximate analytic method, including the Debye - Waller fac- 
tor, in the density range 10'^ g/cva? < p < 10^^ g/cvo? , which corresponds to the cores of 
white dwarfs and outer crusts of neutron stars. Monte Carlo simulations confirmed the 
validity of the results of Itoh et al. (1984, 1993) obtained by the approximate method. 
The analytic approach allowed Baiko and Yakovlev (1995) to fit the kinetic coefficients by 
simple equations, valid for any nuclear composition, and to extend the results to the case 
of the face-centered-cubic Coulomb crystals, which could also exist in dense stellar matter 
but which had not been considered earlier. 

The aim of the present article is to study the kinetic coefficients at higher densities 
10^^ g/cw? < p < 10^^ g/cm^. This range corresponds to the inner crust of a neutron star 
where finite sizes of atomic nuclei are important. Combining the results of the present 
article and the analytic fits of Baiko and Yakovlev (1995) for lower densities, we will obtain 
simple fit expressions valid equally for the outer and inner crusts of neutron stars with 
any nuclear composition. 

2 GENERAL RELATIONSHIPS 

Consider cold matter in the density range 10^1 g/cm^ < p < 10^^ g/cm^. At p < pd 
{pd ~ (4 — 6) X 10^^ g/cm^ is the neutron drip density) matter consists of nuclei {Z,A) 
and strongly degenerate, ultrarelativistic and almost ideal electrons. We assume that, 
at fixed p and T, there are nuclei of single species. Highly energetic electrons induce /? 
captures and, therefore, large neutron excess of atomic nuclei (e.g., Shapiro and Teukolsky, 
1983). At higher densities, Pd < P ^ 10^"^ g/cm^, matter contains also free neutrons in 
addition to the nuclei and electrons. With the grows of p, the fraction of the neutrons 
increases. At p > 10^^ g/cm^ the nuclei may become nonspherical and/or form clusters 
(Lorenz et al., 1993). The nuclear composition depends on pre-history of the neutron star 
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crust. For instance, the composition of ground state matter (Negele and Vautherin, 1973; 
Haensel and Pichon, 1994) differs noticeably from that of accreted matter (Haensel and 
Zdunik, 1990 a, b). 

The properties of matter depend also on temperature T. At low T the nuclei (ions) 
form Coulomb crystal. One commonly studies body-centered-cubic crystals which are 
most tightly bound. However the difference of binding energies of the body-centered-cubic 
and face-centered-cubic crystals is very small, and one cannot exclude the presence of 
face-centered-cubic crystals in the neutron star crust (e.g., Baiko and Yakovlev, 1995). 
The classical body-centered-cubic crystal melts when the ion coupling parameter T = 
Z'^e'^ /{akBT) reaches the critical value T = Tm ~ 172 (Nagara et al., 1987). In this case 
a = [3/(47rni)]^/^ is the mean inter-ion distance, rii is the number density of the nuclei, 
and fee is the Boltzmann constant. Accordingly, the melting temperature is (Figures 1 
and 2) 

^ « 1.323 X lO^zVs ( P^) ^Jl K, (1) 

where /le is the number of baryons per one electron, pQ is density in units of 10^ g/cm^. 

The thermal and electric conduction in crystalline matter (T < Tm) is mainly provided 
by the electrons due to the electron-phonon scattering. Under astrophysical conditions, 
the normal and Umklapp scattering processes are allowed (e.g, Raikh and Yakovlev, 1982). 
The scattering regime is determined by the relationship between the temperature T and 
the ion plasma temperature (Figures 1 and 2) 



Tp = 7.832 xlO%/4^K, (2) 

kB V "^l^e 

where Wp = \J '^'nZ'^e^nij'mi is the ion plasma frequency and is the ion mass. Note 
that the Debye temperature of the crystal is Td = 0.45rp (Carr, 1961). If T > Tp, 
many thermal phonons are excited, and the electron scattering can be treated classically 
(as caused by thermal lattice vibrations). At T ^ Tp the number of thermal phonons 
is strongly reduced, and the scattering is essentially quantum (involves single phonons). 
In the both cases the Umklapp processes dominate the normal ones. However, at very 
low temperatures T ^ Tu (Figures 1 and 2), the Umklapp processes are frozen out, 
and the normal processes become more significant (e.g., Raikh and Yakovlev, 1982). Here, 
Tu ~ TpZ^/^e^/ (Tifp), ^ c being the electron Fermi velocity. We shall restrict ourselves 
to the case of T > Tu, which is most important for applications. In this case, one can use 
the free electron approximation. 

We express the thermal conductivity k and electric conductivity a through the effective 
electron collision frequencies z^^ and v^j: 

, . ''y^"-^4.04xlO-V,5reflg^'| 

om^u^ \ I cm s K 
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^ = ^ « 1.49 X lO^V/? ( 1 -. (3) 

m^Ua y J s 

In this case x = pF/{mcc), fi = v-p/c = x/\/l + x^, m* = meVl + x^, and is the 
electron Fermi momentum. 

The frequencies u,^ and v^y are conveniently expressed (Yakovlev and Urpin, 1980) 
through the dimensionless functions and F^: 

va^n = ^ ^ 0.955 X 10i5^F,,« s-^ (4) 

For calculating these functions, we will use the approximate formulae derived by Baiko 
and Yakovlev (1995). The general formalism is independent of lattice type. The cases of 
the body-centered-cubic and face-centered-cubic crystals differ through numerical factors 
determined by phonon spectrum. According to Baiko and Yakovlev (1995), 

F^ = Go{t)Ko, F^ = F^ + G2{t)(sK2-^Koy (5) 



-2W 



2W 



J Qmin Q 



q |e(g)|2 \^ 4p| J 

Here, is the momentum transferred by an electron due to emission or absorption of 
a phonon, u = [hq / {2pf)]'^ , t = T/Tp is the dimensionless temperature, e~'^^ is the 
Debye — Waller factor, f(q) is the nuclear formfactor, e{q) is the static longitudinal 
dielectric function of the degenerate electron gas (Jancovici, 1962). The lower limit of 
integration over q {q = Qmin) restricts the integration domain where the Umklapp processes 
are operative (e.g., Raikh and Yakovlev, 1982); this limit is set equal to the equivalent 
radius of the Brillouin zone q^am = (GTT^ni)^/^; uq = [%min/(2pF)]^ = 1/(4Z)2/^. 

The functions Go and G2 in Equations (5) are determined by the phonon spectrum. 
They have been calculated and fitted by Yakovlev and Urpin (1980), Raikh and Yakovlev 
(1982), and also by Baiko and Yakovlev (1995): 

Go(t) 



^2(0 - 2/+2 , ^372• (''') 

7r^(i^ + a2) ' 
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In this case U-2, ao and 02 are the constants determined by the lattice type and U-2 is a 
frequency moment of the phonon spectrum. For the body-centered-cubic crystal, we have 
u^^^''^ = 13.00, a^Q"""^ = 0.0174 and of ''''^ = 0.0118. For the face-centered-cubic crystal, 

^^2"^ = 28.80, a'-Q^'''^ = 0.00505 and a^^"''^ = 0.00461. 

As shown by Baiko and Yakovlev (1995), the Debye - Waller parameter can be fitted 
(with the mean error of ^ 1%) as 

2W{q) = au, (8) 
a = ao Qu_ie-9 ioo* + to_2 



X 



where U-i is another frequency moment of the phonon spectrum, u^^i^^ = 2.800, u^l^i^ = 
4.03. 

Since we consider the ultrarelativistic electron gas, it is sufficient to set a; S> 1, /3 = 1. 

For the densities of interest, the sizes of atomic nuclei become comparable to the intcr- 
ion distance a (see above). Contrary to the case of p < 10^^ g/cm^, considered by Baiko 
and Yakovlev (1995), the nuclear formfactor f{q) becomes important in the integrals (6). 
We will not analyze the nonspherical nuclei which may appear at p > 10^^ g/ cm^ (Lorenz 
et al., 1993), but restrict ourselves to a standard model of spherical nuclei with the uniform 
proton core of radius rc- 

3 

/(«) = 7 [sin(grc) - grc cos(grc)] , (9) 



[qr, 



Vr. ... „.l/3 



qr^ = Vu-^ {l8nZ) 



3 NUMERICAL RESULTS AND THEIR APPROXIMA- 
TION 

Calculation of the thermal and electric conductivities reduces to evaluating the integrals 
Kq and K2 from Equations (6). For the conditions of study, the functions Ko and K2 
depend on three parameters: on the quantity a in the Debye — Waller factor (8), on the 
nuclear charge number Z and on the parameter g = rc/a in the nuclear formfactor (9). We 
have calculated the integrals (6) at Z = 20,40,60, a = 0.04, 0.12, 0.4, 1.2, 4, 12, and g = 
0, 0.1, 0.2, 0.3, 0.4. Analyzing the properties of matter at 10^^ g/cm^ ^ P ^ 10^^ g/cm^, 
one can see (Section 4) that the adopted grid covers all possible parameter domain. 

Let us fit the newly calculated values of K^r and K^, basing on the analytic expressions 
for F(j and at p < lO^''' g/cm^ obtained by Baiko and Yakovlev (1995). First of all, note 
that the electric screening is rather insignificant. As shown by Baiko and Yakovlev (1995), 
one can set €{q)=l in the integrals (6) and include the weak screening effect by shifting 



6 



the integration limit uq to ui = (4Z)~^/^ + Ug, where Ug = [?i/ctf/(2pf)]^ = /{'kTivp) k, 
l/(1377r/3) and A;tf is the inverse length of screening of a charge by the degenerate plasma 
electrons. 

First consider the small-size nuclei {g <C 1). Expanding the formfactor up to the terms 
oc g^, we obtain |/(gf)P « 1 - 0.2(grc)^. Then the integrals (6) are taken: 



Here, we have introduced the functions 



$1 



(10) 



Sk-i — P Sk, Sk 



1 k —cxu 

auu e 



(11) 



Ul 

1 



5_i = E{w)-E{a), So = - (e-'" - e-") , 
Si = ^ [e-^(w; + l)-e-"(a + l)] , 



So 



(12) 



where w = aui, E{x) is the integral exponent, which can be easily calculated using, for 
instance, Equation (21) of the article of Baiko and Yakovlev (1995). 

The natural generalization of Equation (10) to the case of large-size nuclei can be 
provided by the following expressions: 



^0 
K2 



2#i 



2 " 



1 + 



(187rZ)2/3ff2^2" 

5$iPo 
5$o-P2 



-P-2 



(13) 
(14) 



where Pq and P2 are the fit parameters. These expressions describe quite accurately our 
numerical results with Pq = 4.787 - 0.0346 Z and P2 = 2.729 - 0.0204 The mean fit 
error of Equation (13) is 2%; the maximum error 4.7% takes place at a = 0.04, Z = QO 
and g = 0.3. The mean fit error of (14) is 2.3%, and the maximum error 10% takes place 
at a = 12, Z = 60 and g = 0.4. 

If 3 <C 1, Equations (13) and (14) reproduce the fit expressions (18) - (27) of the article 
of Baiko and Yakovlev (1995). The latter equations are valid at rather low densities 10^ 
g/cm^ < p < 10^^ g/cm^, at which the finite sizes of the nuclei are unimportant. Thus 
Equations (13) and (14) cover very large density range 10^ g/cm^ ^ P ^ lO^'* g/cm^ (outer 
and inner neutron star crusts) and are valid for any nuclear composition. 
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4 ELECTRIC AND THERMAL CONDUCTIVITIES 

Let us analyze the electric and thermal conductivities of electrons in crystalline matter of 
a neutron star crust. The nuclear composition of matter is not definitely known. For illus- 
tration, consider two well established models: ground state matter and accreted matter. 

Ground state matter corresponds to the minimum of free energy per one baryon. This 
matter can be produced from hot dense matter of a newly born neutron star in the course 
of the subsequent stellar cooling. At densities p < pd ^ 4 x 10^"*^ g/cm'^ (below the 
neutron drip, Section 2) we will use the composition of the ground state matter calculated 
by Haensel and Pichon (1994) on the basis of new laboratory measurements of masses of 
nuclei with large neutron excess. At p > we will use the results of Negele and Vautherin 
(1973) derived with the aid of a modified Hartree-Fock method. 

The nuclear composition of accreted matter is determined (Haensel and Zdunik, 1990a, 
b) by nuclear transmutations (pycnonuclear reactions, /^-captures, emission and absorption 
of neutrons, etc) in cold dense accreted matter sinking within the star. The results are 
practically insensitive to the accretion rate and initial chemical composition. The accreted 
matter is assumed to burn into ^^Fe at p < 10^ g/cm^. 

Both models of matter are formally calculated at T = 0, but they are actually valid at 
r < 4 X 10^ K (as long as the thermal effects do not influence the properties of nuclei). It 
is adopted that single species of nuclei is present at every pressure (density). This leads 
to jumps of the nuclear composition (A, Z) with the growth of pressure (density) and to 
jumps of typical temperatures depicted in Figures 1 and 2. At p < 10^ g/cm^ matter 
consists of ^^Fe in both models. At higher densities, the compositions of ground state and 
accreted matters are noticeably different. The accreted matter consists of lighter nuclei 
with lower Z. For instance, at p = 4 x 10^^ g/cm^ the accreted matter contains neutron 
rich magnesium nuclei, A = 4A, Z = 12, while the ground state matter contains tin 
nuclei, A = 159. Z = 50. The difference of compositions affects the thermodynamic and 
kinetic properties of matter. For example, at 10^^ g/cm^ ^ P ^ 10^^ g/cm^ the melting 
temperature of the accreted matter is a factor of 3 — 10 lower, than that of the ground 
state matter. 

One should know the proton core radius Vc, for calculating the thermal and electric 
conductivities. Following Itoh et al. (1984) we set: rc = 1.15 ^^^^ fm at p < pd', Vc = 
1.83 fm at p > pd- It is easy to verify that g = Vc/a < 0.2 for both models. Figure 
3 shows the functions F^j and versus T/Tp for a body-centered-cubic lattice of nuclei 
A = 159, Z = 50 (p = 4 X 10^^ g/cm^). The curves are qualitatively the same as in 
matter of lower density (p < 10"*^^ g/cm'^; see Baiko and Yakovlev, 1995). The allowance 
for the nuclear size decreases F^-^k and the effective electron collision frequencies Va,K-, i-e-, 
increases the electric and thermal conductivities. However, the effect is not strong, about 
(20 — 30)%. This is clearly seen from the integrals (6): the nuclear fornifactor reduces the 
integrands with the growth of g only at large dimensionless momentum transfers m ~ 1 
(i.e., at large-angle scattering), at which the integrands themselves are sufficiently small. 
The latter smallness is due to the Debye-Waller factor (which suppresses collisions with 
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large momentum transfers) and due to the factor (1 — P^u), that describes reduction of 
backward scattering of relativistic electrons. 

Figure 4 shows the dependence of F^^^ on T/Tp for the same conditions as in Figure 3, 
but for the face-centered-cubic lattice. As in a lower-density matter (Baiko and Yakovlev, 
1995), the lattice type affects strongly the kinetic coefficients. The finite sizes of the nuclei 
again influence the result weakly, as in the body-centered-cubic lattice. 

Figures 5 and 6 display the density dependence of the electric and thermal conductiv- 
ities for the ground state and accreted matters at T = 10^ and 5 x 10^ K. At low p the 
curves are broken in melting points (Figures 1 and 2). The discontinuities of the kinetic 
coefficients are associated with the jumps of nuclear compositions (see above). At 10^^ 
g/cm^ < p < 10^^ g/cm^ the electric and thermal conductivities of the accreted matter 
are a factor of 2 — 3 higher, than those of the ground state matter. This is naturally 
explained by lower nucleus charge (and lower electron collision frequencies) in the accreted 
matter. 

Note that earlier the thermal and electric conductivities of the ground state matter 
in the inner crust of a neutron star was calculated and fitted by Itoh et al. (1984). Our 
results agree mainly with the results of these authors. However, there is a 30 — 50 % 
difference in the thermal conductivity at densities p > 10^^ g/cm^. Let us emphasize that 
our fit expressions are much simpler than those of Itoh et al. (1984), and they are valid 
for matter with any nuclear composition. 



5 THERMAL DRIFT OF THE MAGNETIC FIELD 



If the magnetic field is present in the stellar plasma, various thermal magnetic effects may 
operate (Urpin and Yakovlev, 1980). For instance, the Hall component of the thermopower 
affects the thermal flux, that emerges from the star, and induces the thermal drift of the 
magnetic field. If the field is not too large, so that the electron gas is weakly magnetized 
{lob t <^ 1, where u>b = eB/{cm*)), one can neglect the back reaction of the field onto 
the thermal flux. According to Equations (5) and (7), at T > Tp one has Ffj ~ F^. This 
allows us to introduce the electron relaxation time ro(e) = = i^^^, which is the same 
for the thermal and electric conductions {s being the electron energy). Under formulated 
conditions, the drift velocity u is parallel to the thermal flux Q (Urpin and Yakovlev, 1980) 



u 




(15) 



1 dln{To{s)e-^) 

2 dlnp 



P=Pf 



(16) 



Here, we have assumed that Q = cqT^, where Tg is the effective temperature of the stellar 
surface and ctq is the Stefan - Boltzmann constant. The factor 77 ~ 1 and can change sign. 
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that corresponds to different thermal drift directions. Taking the derivative, we have 



(1-/3V 



I/(2pfA)P 
|e(2pFM)|' 



+ r du 

Juo 



\ml 



ue 



■ (17) 



The dashed hne in Figures 1 and 2 shows those densities and temperatures at which rj = 
0. One has 77 > above this Hne, i.e., the thermal drift is directed outward, and rj < 
below the line, so that the magnetic field drifts inward. 

Note that Equations (15) and (16) are actually valid only at Tp < T < Tm, and the 
dashed line is broken when T reaches with the growth of p. It is seen that the thermal 
drift can be directed either inward or outward in the above temperature range and at 
densities up to lO^'^ - 10^^ g/cm^. At higher densities and T > Tp, the factor 77 remains 
positive for the ground state matter. As for the accreted matter, Tm exceeds Tp at these 
densities, and Equation (16) becomes invalid. Calculation of the thermal drift at T < Tp 
is outside the scope of the present article. 



6 CONCLUSIONS 

We have calculated the thermal and electric conductivities of relativistic degenerate elec- 
trons due to electron-phonon scattering in Coulomb crystals of atomic nuclei at densities 
10^^ g/cm^ < p < 10^^ g/cm^, corresponding to the inner crust of a neutron star. Calcu- 
lations are done with the aid of the approximate analytic method (Baiko and Yakovlev, 
1995) taking exact spectrum of phonons and the Debye-Waller factor into account. In 
addition, wc have taken into consideration finite sizes of atomic nuclei. 

The results are fitted by a few simple equations which reproduce also the fit expressions 
obtained by Baiko and Yakovlev (1995) for lower densities. Thus we have obtained the 
unified fits which describe the thermal and electric conductivities in the wide density range 
10^ g/cm^ ^ p < 10^^ g/cm^ (in the inner and outer crusts of neutron stars) below the 
melting temperature but above the temperature where the Umklapp processes are frozen 
out, for the body-centered-cubic and face-centered-cubic crystals. Note that the same 
equations are valid in the presence of a weak magnetic field {lob t ^ 1). 

Our results can be useful for numerical simulations of thermal evolution of neutron 
stars (cooling, nuclear burning of accreted matter) and evolution of their magnetic fields 
(ohmic dissipation, generation due to thermal magnetic effects). 
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Figure 1: Density — temperature diagram for the ground state matter: Tm is the melting 
temperature, Tp is the ion plasma temperature, Tu is the freezing temperature of Umklapp 
processes. Dashed line corresponds to the conditions at which the thermal drift velocity 
equals zero. 
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Figure 2: Same as in Figure 1 but for the accreted matter. 
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Figure 3: Factors F„ (curves 1) and (curves 2) versus t = T/T-p for the body-centered- 

cubic and face-centered-cubic crystals composed of ^"''^Sn nuclei at density p = A x 10^^ 
g/cm^ with account for the finite nucleus size {g = 0.15, solid lines) and for the point-like 
nuclei {g = 0, dashes). 
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Figure 4: Same as in Figure 3 but for the face-centered-cubic crystal. 
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Figure 5: Electric conductivity versus density for the ground state (solid lines) and 
accreted (dashes) matters at T = 10^ K (curves 1) and T = 5 x 10^ K (curves 2). 
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Figure 6: Thermal conductivity versus density for the same conditions as in Figure 5. 



